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1. Introduction 

Most applications of diffraction theory in optics are 
based on the scalar approximation of Huygens' princi- 
ple, and it is usually understood that polarization is 
ignored. A notable exception to this practice is the clas- 
sical pseudo-vectorial theory formulated by Rayleigh 
[1,2] and Bouwkamp [3] for the diffraction of plane- 
polarized light by small metallic apertures as illustrated 
in Fig. 1. The results obtained by these authors implied 
large polarization effects in the case of aperture sizes 
smaller than the wavelength of light, X. Specifically 
Rayleigh and Bouwkamp predicted that unpolarized 
light diffracted by very narrow aperture will become 
strongly ^-polarized, as indicated by the ratio T^/rj^ 
of the transmission coefficients of circular apertures 
and slits illuminated by normally incident^- or s-polar- 
ized light plotted in Fig. 2 for the range < kw < 0.75, 
where 2w is the aperture width and k = 2k/ X is the cir- 
cular wave number. Regarding the case of slits, 



Rayleigh commented: "It appeared that if the width of 
the slit is very small in comparison with the wave- 
length, there is a much more free passage when the 
electric vector is perpendicular to the slit than when it 
is parallel to the slit, so that unpolarized light incident 
upon the screen will, after passage, appear polarized in 
the former manner." While this may be plausible as far 
as slits are concerned, it seems impossible that a circu- 
lar aperture can behave in the same manner because 
normal incidence is assumed so that symmetry would 
dictate a complete absence of polarization. Apparently 
Rayleigh had similar doubts. He remarked that the p- 
solution has "no simple application" in the case of cir- 
cular apertures, but nonetheless the discrepancy is dis- 
turbing as the Rayleigh-Bouwkamp theory is common- 
ly referred to as a "rigorous" theory, and in this connec- 
tion it should also be noted that later authors, such as 
Levine and Schwinger [4,5], considered only the s- 
solution and, in doing so, made no reference to polar- 
ization. In the present paper the Rayleigh-Bouwkamp 
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Fig. 1. The Rayleigh-Bouwkamp diffraction problem. A very narrow aperture Q$f in a metal- 
lic screen d? 9 which is located in the xy-plane is illuminated by a monochromatic, p- or ^-polar- 
ized plane wave e 1 z , and solutions are constructed so that the optical field is everywhere con- 
tinuous with continuous first derivatives. 
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Fig. 2. Polarization ratios, p = %/%» of circular apertures and slits vs. aperture size kw according to 
Bouwkamp [3, p.71ff]. 
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theory will be reexamined with respect to its physical 
significance, mathematical rigor, and accuracy 

Rayleigh considered a monochromatic plane wave 
of unit irradiance, u = e lfe , which is incident from the 
half space z < on a circular aperture or slit osf con- 
tained in an infinite, infinitesimally thin, opaque screen 
d^that lies in the xy-plane of a cartesian coordinate sys- 
tem (Fig. 1). He also assumed that the diffracted field 
in the half space z > obeys the Rayleigh-Sommerfeld 
integral equations, 



„«(P) = -±fdQ*^®^ z,0, (1 
?- J V dz QP V 



2n 



a) 



« W (P) = 



J_ 

2% 



JdQ« w (Q) 



a (e kQP ^ 



dz 



QP 



z>0, (lb) 



where Q = (£, 7], 0) is a point in o2^ 9 P = (x>y> z) is an 
arbitrary point in space, QP is the distance between 
them, u^ ,s) are scalar wave functions, and d u^\Q)/d z 
and u {s \Q) are zero on the opaque portion of <§? as if the 
screen is a perfect metallic reflector and the incident 
light is p- or s-polarized. Furthermore, Rayleigh 
assumed that field on the source side of <$* is a modi- 
fied geometrical field of the form 

u (p ' s) (P) = c ikz ±^ +u ( _ p ' s) (P), z<0, (lc) 

and hence the objective of the theory is to find the 
unknown quantities u^XP) and u^XP) so that the 
overall field is continuously differentiable in the aper- 
ture plane. 

From Eq. (1 c) it follows that d u^idz and u (s) will be 
continuous for z = if 

^WO) .MM, M , W , ) =B « ( ^o), 



dz 



dz 



(2a) 



so that 

u["' s) (x,y,z) = + u ip ' s) (x,y,-z), z>0, (2b) 

u ip ' s) (x,y,z)=c ik ±c- ih +u ip ' s) (x,y,-z), z<0. (2c) 

Thus, « w = 2 - « w , 3« w /3z = 2i£ - a« w /az when z = 0, 
so that w (p) and du (s) /dz will also be continuous if 



ik dz 



(2d) 



The conditions [Eq. (2d)] represent the underlying 
integral equations of the Rayleigh-Bouwkamp and 
other so-called "rigorous" diffraction theories which 
are intended to solve the wave equation in the presence 
of prescribed boundary conditions. To this author's 
knowledge, all existing solutions with the sole excep- 
tion of Sommerfeld's half-plane theory are approxima- 
tions of some sort and in the cases considered here 
these approximations take the form of ascending power 
series in kw which were derived on the assumption that 
the unknown aperture distributions in Eqs. (la,b) can 
be expressed in the form 



1 du ip) (Q) 
ik dz 



l 

kw 



^ + {kwfpW + ... 
W l 



u {s) (Q) = -ikw [s x W + (kwys 2 W 5 +...], 



(3) 



where W = y]l- (ql w) 2 , q is the distance of Q from the 
coordinate origin, and the coefficients p n and s n are to 
be chosen so that Eqs. (2d) are satisfied. The determi- 
nation of these coefficients is tedious and will not be 
repeated here as the results have been detailed in Ref. 
[3]. In the following the corresponding solutions will 
be denoted by the subscript 'RB' and it will be suffi- 
cient to consider only the initial terms of the expansions 
[Eq. (3)] by assuming, as Rayleigh did originally, that 
kw -> 0. 



2. Circular Apertures 

For circular apertures it is appropriate to use polar 
coordinates inside the aperture and in the plane of 
observation, so that 

P = (p cos<p,p sin<p, z), p 2 =x 2 +y 2 , (4a) 

Q = (q cosx,q sin^, 0), q 1 = «f +7] 2 , (4b) 

dQ = qdqdx, QP 2 =z 2 +p 2 +q 2 -2pqcos(<p-x). (4c) 

The initial terms of the trial solutions of Eq. (3) were 
determined by Rayleigh and Bouwkamp as p = s l = 
-21k and thus the first-order solutions for circular aper- 
tures are 



1 w 2S ne- ikQP 

n o o y]w-q 2 t 



q 2 QP 



- lcw->0, z>0, 



(4d) 
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2.1 Fresnel Approximation 

The evaluation of Eqs. (4d,e) in the Fresnel limit is 
simple. When z » X, and thus z » q, the dependence 
of QP on the aperture coordinates can be ignored, so 
that 



QP~*Jz 2 +p 2 = r, - = cos0 



and 
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where q/w = sina was substituted in both integrals. 
Hence, the total radiant fluxes transmitted by the aper- 
ture and incident on the inside of a large hemisphere 
with radius r are given by 

O™ = 2w 2 *[ d<9sin<9| <> | 2 = — , ( 5( j) 

o * 

jc/2 ol4 6 n/2 

<B w =2jcr 2 f d0sin0|w£>| 2 =-5^- fd0sin0cos 2 



9rc 



8Tw_ 
27rc 



(5e) 



and the corresponding transmission coefficients of the 
aperture are 

_0^___8_ w _ <B (J) _ 8(frv) 4 






o ^ 



-, Aw^O, (5f) 



where O = 7iw 2 is the geometrical flux incident on the 
aperture. 

These results can now be checked for consistency 
with Eqs. (2d) and (3) by using an alternative expres- 
sion [3,6] for transmission coefficients, 



r' = UdQ\u(Q)\ 



(6a) 



where A is the area of the aperture and w(Q) is the field 
inside it. According to Eq. (2d) we expect w^(Q) = 1 
and therefore 



whereas from Eq. (3) and s 1 = -21% it follows that 

\2 2% 1 xi/7_..\2 



(6b) 



TT\A? J J ft 

(6c) 
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The ^-coefficients of Eqs. (5f) and (6b) agree reason- 
ably well. On the other hand, there is a considerable 
difference between the corresponding ^-coefficients. 

2.2 Near-Field Solution 

In the immediate proximity of the aperture plane the 
exponential e lkQP is effectively constant and equal to 1 
inside the aperture, and Eqs, (4d,e) can be simplified as 
follows, 



1 } qdq 2 ?d X 






-*-, z«X (7a) 



(7b) 



As these wave functions must be rotationally sym- 
metrical about the z-axis it suffices to evaluate the %- 
integrals appearing in these expressions for <p = 0. 
Substituting a={%- rc)/2, we find 

QP 2 = z 2 +p 2 +q 2 -2pqcosx=—0 ~ msin 2 a), (8a) 

m 



m = 



*pq 



z 2 +(p + q) 2 



(8b) 



and hence 

2 



[ — = 2 J — | — — = 2 A — K(m), (8c) 

{QP "\pq{{\-msm 2 af 2 Vp<7 

2 f dj m ^r da _ tc/w 
o Q p2 P<7 o ] " m ™ 2 a 2pq^T^} 



*, (8d) 



Volume 111, Number 1 , January-February 2006 

Journal of Research of the National Institute of Standards and Technology 



\ 3/2 71/2 

H-) f- 



da 



QP 3 l\pq ) { (l-msin 2 a) 3/2 



1 f m 



\3/2 



E(iw) 



- m 
(8e) 



where K(m) and E(m) are the complete elliptical inte- 
grals of the first and second kinds [7,8]. 

Upon substitution of Eq. (8c), Eq. (7a) is reduced to 
the single integral 



<(P): 



\ d(q I w)Jmq i wK(m) 

j , , z«A, 



n 2 yjp/w{ sjl-iq/w) 2 



(9a) 



where relative coordinates are used. This expression 
was evaluated by numerical integration to analyze its 
behavior in the aperture plane, using a polynomial 
approximation for K(m) given in Ref. [7] and assuming 
z/w = as well as z/w = 0.01 in order to find du^/dz by 
numerical differentiation. The results of these computa- 
tions are listed in the left-hand portion of Table 1, 
which shows that they are in satisfactory agreement 
with the values expected from Eqs. (2d) and (4d). 

The corresponding expressions for u^ are obtained 
from Eqs. (7b) and (8d) and can be expressed as fol- 
lows, 



«S(P) = fe 
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where 
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and E(m) can also be expressed by a polynomial 
approximation [7]. It was ascertained that the integrals 
A and B are everywhere finite, and hence it follows 
from the first Eq. (10a) that the aperture values of u^ 
are identically equal to zero. This is in disagreement 
with Eq. (4e), and it was also found that the computed 
aperture values of (l/ik)du^/dz do not satisfy the sec- 
ond integral equation, [Eq. (2d)]. These discrepancies 
are documented in the right-hand portion of Table 1 . 



3. Slits 

The second case considered by Rayleigh and 
Bouwkamp is that of an infinitely long, narrow slit of 
width 2tv which is centered on the j-axis of a rectangu- 
lar coordinate system. As the resulting diffraction pat- 
tern will be independent of y, it is sufficient to compute 
its variation in the xz-plane by assuming 



P = (jc,0,z), Q= (&*7,0), 
and then Eqs. (la,b) are reduced to 



(lla) 



<(P)="T L Jd^ 



Jk^-x) 2 +Tl 2 +z 2 



1 7..ftfl(g,0)j d e: 
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«S(P)=~Jd5«e(5,o)^^ 
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^ -TV P 



where the infinite ^-integral in Eq. (10b) was evaluat- 
ed in terms of the Hankel function H^\p) [8] and the 
relationship dH^/d^ = -H^/J) was used in Eq. 
(10c). According to Rayleigh and Bouwkamp the initial 
terms of the series expansions [Eq. (3)] arejt? = l/(o - 
in/2) and s 1 = 1 where and o = h\(ykw/4) and y= 0.577.. 



Table 1. Aperture values of wave functions and normal derivatives for circular apertures of width kw = 0.1, as computed from Eqs. (9a) and 
(10a,b,c) and expected from Eqs. (2d) and (3) 



p/w 



«Scs) 



computed 



expected 



(l/ik)du^Qydz 



computed 



expected 



computed expected 



(l/ik)du ( ^(Q)/dz 



computed 



expected 



0.01 


0.99 


L -0.62 


-0.66 





0.064 i 


0.5 i 1 


0.25 


0.99 


L -0.64 


-0.66 





0.062 i 


-30 + 0.5 i 1 


0.50 


0.99 


L -0.72 


-0.74 





0.055 i 


-216 + 0.41 1 


0.75 


0.99 


L -0.94 


-0.96 





0.042 i 


-558 + 0.3 i 1 


0.99 


0.98 


L -3.91 


-4.51 





0.009 i 


-276 1 
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is Euler's constant, and hence the first-order solutions 
for slits are 

«4g(P) = - H *&< )(P) , ^->o, z>o, (iid) 
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3.1 Fresnel Approximation 



At large distances from the slit (z » w), the phase 
difference {} defined in Eq. (l ib) is ef fectively a large 
constant and equal to kr = Wx 2 + z 2 , and the Hankel 
functions in Eqs. (lld,e) can be replaced by their 
asymptotic expansions 



iTO)' 



2* \(1rr-TrlA~\ ^Q \P ) 



,i(Jfa--7i/4) 



n(kry 



k i(Jfa--37i/4) 
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so that 
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o 
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As in Sec. 2. 1 alternative expressions for these trans- 
mission coefficients can be obtained from Eq. (6a), 
which in this case is equivalent to 



^ TV 

T' = -fd£|«(4)| 2 - 

vv 



(13a) 



Hence we find, using Eqs. (2d) and (3) with s 1 = 1, 

(0' = 1> (13b) 

(r^y = kw\^[\-aiwf]=^^, kw^O. (13c) 



Both of these results disagree with the corresponding 
values in Eq. (12e). 

3.2 Near-Field Solution 

The integrals [Eq. (lld,e)] can be evaluated without 
simplifying assumptions for arbitrary values of z, and it 
was found that they are finite in the aperture plane in 
spite of the singularities of H< 1} 03) and H^/J)//? for 
0=0. 

For ^-polarized light we find 

_ 1 l , dg/w){[7 I J (ffl /2+oY (^)]-i[oJ (^)-7 I Y (^)/2]} 

Oil/ J 



2M* 



VH£/w) 2 



(14) 



where cos0 = z/r and the two integrals were again eval- 
uated by substituting !~/w = sina. These expressions can 
now be used to compute the quantities 



III z 

AO^ s) = rAy J d<9 \ u 



(p,s) .2 

■SB I s 



which represent the radiant flux emerging from a slit 
segment of length Ay and falling on the corresponding 
annular segment of a semicylinder of radius r which is 
centered on the^-axis. As the geometrical flux incident 
on the slit segment is A<3? = 2wAy, this leads at once to 
the following expressions for the transmission coeffi- 
cients of the slit, 



where J n and Y n are Bessel functions of the first and 
second kinds andM= o 2 + n 2 /4. The integral [Eq. (14)] 
was evaluated numerically for kw = 0. 1 and z = as 
well as z = 0.01 w, yielding the results for u^ and 
(12d) du^idz shown in the left-hand portion of Table 2. The 

table also shows the corresponding values expected 
from Eqs. (2d) and (3), and the agreement was judged 
satisfactory in view of the approximations made in 
Rayleigh's derivations. 

The corresponding expression obtained from Eq. 
(lie) is 



«S(P) = - 



k 2 wz 



(C + iZ)), 



(15a) 
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Table 2. Aperture values of wave functions and normal derivatives for slits of width kw : 
expected from Eqs. (2d) and (3) 



0.1, as computed from Eqs. (14) and (15a,b,c) and 



plw 



computed 



«rbOS) 



expected 



(i/i*)aw§key9z 

computed expected 



» 



(Q) 



computed 



expected 



(Uik)du { ^(Q)/dz 



computed 



expected 



0.01 


0.76 


L 0.75 + 2.0 i 


0.76 + 2.11 





-0.10 i 


8.6-381 1 


0.25 


0.76 


L 0.78 + 2.1 i 


0.79 + 2.1 i 





-0.10 i 


8.3-37 1 1 


0.50 


0.76 


1 0.87 + 2.3 i 


0.89 + 2.4 i 





-0.09 i 


7.4-33 1 1 


0.75 


0.76 


L 1.1+3.H 


1.2 + 3.1 i 





-0.05 i 


5.8-25 i 1 


0.99 


0.75 


L 4.7+131 


5.5 + 15 i 





-0.01 i 


1.6-4 i 1 



where 



_\d^/w)^-^/w) 2 J l (p) 



C=J 



I d(<g/w)Vl-(<g/w) 2 ¥,(/?) ^ 



It was found that the integrals [Eq. (15b)] are finite in 
the aperture plane, so that according to Eq. (15a) the 
aperture values of u$ are identically zero and thus do 
not agree with the assumed values of Eq. (3). Equation 
(15a) also shows that 



1 



_3#) 
\k dz 



=f(/>-iQ 



if z=0, 



(15c) 



which obviously disagrees with Eq. (2d). The various 
numerical results obtained for u$ and du^/dz are list- 
ed on the right-hand side of Table 2. 



4. Conclusion 

Based on the above-mentioned results and illustrated 
by the numerical examples shown in Tables 1 and 2 
there can be no doubt that the Rayleigh-Bouwkamp 
theory is mathematically flawed and fails to give a 
credible account of the diffraction of polarized light. 
The reasons for this conclusion are twofold. 

The near-field solutions derived in Sees 2.2 and 3.2 
for ^-polarized incident light (dufiz = on d? 9 ) are in 
approximate agreement with the conditions in Eqs. (2d) 
and (3) that were assumed in their derivation, and thus 
appear to be free of contradictions. On the other hand, 
the corresponding solutions for ^-polarization (w = 
on d? 9 ) do not obey these conditions at all. This discrep- 
ancy was found for circular apertures and slits alike, 
and thus points to a general problem associated with the 
second integral equation [Eq. (2d)] as such. Bouwkamp 
[3] mentioned that this integral equation [Eq. (2d)] is 
notoriously difficult as it requires differentiation with 



respect to z under the integral sign and can lead to 
divergent solutions. He pointed out that these divergen- 
cies can be eliminated, but this was not attempted in 
this work. In any case, the polarization effects shown in 
Fig. 1 must be dismissed as wrong. 

A second, and most likely uncorrectable, problem 
arises from the fact that the 'polarization effects' pre- 
dicted by the Rayleigh-Bouwkamp theory run counter 
to common experience. As mentioned in the 
Introduction, a polarized field behind a circular aper- 
ture illuminated by normally incident, unpolarized light 
is impossible for reasons of symmetry. Yet the exis- 
tence of such a field is implied by the upper curve in 
Fig. 1 , and the discrepancy persists even when the sus- 
pect transmission coefficient t^ } in Eq. (5f) is replaced 
by the improved coefficient (r^)' in Eq. (6c). 
Additionally, the different values obtained in Sees. 2.1 
and 3.1 for the far- and near-field transmission coeffi- 
cients seemed paradoxical because the polarization of 
light cannot change during free-space propagation. As 
this discrepancy was encountered for both states of 
polarization, it also cannot be blamed on the mathemat- 
ical problem mentioned in the previous paragraph. 

Hence it is concluded that the Rayleigh-Bouwkamp 
theory cannot be relied upon for describing the diffrac- 
tion of polarized light by small apertures. Virtually the 
same determination was made in an earlier paper [6] on 
the feasibility of using the Rayleigh-Sommerfeld 
boundary-value integrals for analyzing diffraction by 
apertures of arbitrary size, and thus the overall conclu- 
sion is that the inability to account for polarization 
effects is an inherent property of scalar diffraction the- 
ories, no matter what boundary conditions are assumed. 
The fact of the matter appears to be that, in scalar dif- 
fraction theory, the boundary conditions du Q /dz = or 
u Q = are necessary and sufficient to satisfy Helmholtz' 
theorem. As it happens, the same conditions also per- 
tain to the metallic reflection of p- or ^-polarized light 
but it does not follow that they are sufficient to ensure 
conformance with electromagnetic theory. 
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In Ref. [6] it was also shown that a modified theory 
of unpolarized diffraction can be formulated in which 
the Rayleigh-Sommerfeld integrals are used to describe 
a continuous, bidirectional flow of energy in the near 
zones on either side of an aperture. The transmission 
coefficients computed in this manner are similar to 
those published by Levine and Schwinger [4,5] for sub- 
wavelength aperture sizes, so that a special theory for 
this region may not be needed. 
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